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Abstract 

The electroweak phase transition and the sphaleron decoupling condition in the minimal super- 
symmetric standard model are revisited taking the latest experimental data into account. The 
light Higgs boson scenario and the ordinary decoupling limit which are classified by the relative 
size between the CP-odd Higgs boson mass and Z boson mass are considered within the context 
of electroweak baryogenesis. We investigate v/T at not only the critical temperature at which the 
effective potential has two degenerate minima but also the nucleation temperature of the critical 
bubbles, where v is a vacuum expectation value of the Higgs boson and T denotes a temperature. 
It is found that v/T at the nucleation temperature can be enhanced by about 10% compared to 
that at the critical temperature. 

We also evaluate the sphaleron decoupling condition including the zero mode factors of the fluc- 
tuations around sphaleron. It is observed that the sphaleron decoupling condition at the nucleation 
temperature is given by v/T > 1.38 for the typical parameter sets. In any phenomenologically al- 
lowed region, v/T at both the critical and nucleation temperatures cannot be large enough to 
satisfy such a sphaleron decoupling condition. 
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I. INTRODUCTION 



The origin of the baryon asymmetry of the Universe (BAU) has been a great puzzle 
in particle physics and cosmology. From the latest observations of the cosmic microwave 
background radiations, the supernovae and the large scale structures, the baryon-to-photon 
ratio is given by [l| 

n B 



— = (4.7 - 6.5) x 1(T 10 (95% C.L.). 



ru 



10|, the singlet-extended MSSMs 



111, 



(1) 

The nonzero BAU can be generated if the following three conditions are fulfilled [2]: (I) 
baryon number (B) violation, (II) C and CP violation, and (III) departure from thermal 
equilibrium. Electroweak baryogenesis is a scenario in which the BAU is explained based only 

n 

on electroweak physics |3j. The standard model (SM) can in principle accommodate above 
three conditions. However, it turns out that the SM fails to generate the BAU due to the 
lack of sufficient CP violation M , and the electroweak phase transition (EWPT) is a smooth 
cross over for > 73 GeV [5j, which makes the condition (III) infeasible. This failure 
motivates us to go beyond the SM. So far, many attempts of electroweak baryogenesis have 
been made in the various new physics models such as the minimal supersymmetric standard 
model (MSSM) 0,33, the two-Higgs-doublet model 
and recently the fourth-generation models also have been paid much attention [12j. Among 
them, the MSSM is a theoretically well-motivated model and a good candidate for successful 
electroweak baryogenesis. 

As far as the strong first order EWPT which is necessary for condition (III) is concerned, 
the light Higgs boson is generically favored. On the other hand, the negative results of 
the Hig gs b oson (s) searches at LEP put the constraints on the Higgs boson masses in the 
MSSM 113 



[yj]. The analysis has been done based on several benchmarks rather than 
scanning whole parameter space. Recently, it is pointed out that the LEP constraints can 
be relaxed in the more general case. In Ref. 15] . it is found that the lightest Higgs boson 
mass can be smaller than the Z boson mass. Such a light Higgs boson can be viable since 
its coupling to Z boson is modified by the additional Higgs bosons and suppressed enough 
to avoid the experimental observed limits. We call a light Higgs boson scenario (LHS) when 
the mass of the lightest Higgs boson is less than 114.4 GeV. This scenario can accommodate 



the 2.3cr excess around 98 GeV Higgs boson mass in the LEP experiments 13j, 1 14] and 



have parameter space that is consistent with the observed dark matter relic density 3] • In 



Ref. [171 ]. it is claimed that the LHS can provide a solution for a little hierarchy problem 
as well. As for the electroweak baryogenesis context, since the neutral Higgs bosons are 
relatively light, they may give a new window to a viable baryogenesis scenario. 

In the electroweak baryogenesis mechanism, a sphaleron process must be decoupled after 
the EWPT in order to avoid the washout of the generated BAU. Because of this requirement, 
the ratio of the Higgs vacuum expectation value (VEV) to the temperature in question should 
be greater than some value. To search for the region where the strong first order EWPT is 
realized, and which is also experimentally allowed, not only the theoretical calculations of 
the zero temperature observables but also those at the finite temperature are mandatory. 
The accuracy of the former calculations are already beyond at least the tree level. On the 
other hand, the sphaleron decoupling condition, which is derived using the tree-level Higgs 
potential at zero temperature, has been frequently used to discuss the strength of the strong 
first order EWPT in the literature. Namely, Vc > Tq is adopted as a practical criterion, 
where Tq is the critical temperature at which the Higgs potential has two degenerate minima 
and vc is the Higgs VEV at To- To obtain more accurate theoretical predictions, the refined 
calculation of the sphaleron decoupling condition is indispensable. The uncertainties of 
the condition is expected to be reduced up to about 10% level if we make the following 
improvements, i.e., the sphaleron solution at finite temperature is used and the contributions 
from the zero modes of the fluctuations around sphaleron at finite temperature are taken 
into account. In addition, the sphaleron decoupling condition must be imposed at the 
temperature Tg lower than Tc, at which the EWPT terminates. In practice, since it is 
difficult to determine this temperature, we substitute T/v, a nucleation temperature of the 
critical bubble, in place of Tg. The critical bubble is defined as the bubble whose surface 
energy and volume energy becomes balanced. Only such bubbles can nucleate and expand 
in the symmetric phase. 

In this article, the possibility of the strong first order EWPT in the LHS and the con- 
ventional decoupling limit is examined with a particular emphasis on the refinement of the 
sphaleron decoupling condition. The analysis is conducted taking the latest experimental 
data into account as well. It is well known that a light stop whose mass is smaller than that 
of top quark is required for the strong first order EWPT. We search for the parameter space 
that is consistent with such a light stop. In the LHS, since the masses of all the Higgs bosons 
are relatively light, typically, < 140 GeV, the B physics observables data can restrict the 
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allowed region. Here, we consider the constraints from the following processes: B u — > tu, 
B — > Xg'-f and B s — > yU + /i~. The first two modes are relevant when the charged Higgs bosons 
are light. The last one, which has not been observed yet, can give a significant limit when 
the neutral Higgs boson is light. 

After finding the phenomenologically allowed region, we investigate the strength of the 
EWPT at not only Tq but also T^. By using the one- loop effective potential at zero and finite 
temperatures, we calculate the sphaleron energy and the contributions from zero modes of 
the fluctuations around sphaleron to evaluate the sphaleron decoupling condition. It is then 
demonstrated whether or not such an improved sphaleron decoupling condition is satisfied 
at Tc and T^. 

This article is organized as follows: In Sec. HIl we present the Lagrangian of the gauge- 
Higgs system for describing the statics and dynamics of the EWPT. The LHS is discussed in 
Sec. IIHI The qualitative evaluation of the EWPT is outlined in Sec. [TV] The critical bubbles 
and the sphaleron decoupling condition are analyzed in Sees. |V] and IVTl respectively. The 
experimental constraints are considered in Sec. IVII1 In Sec. IVIIIl we show the numerical 
results. Section [IX] is devoted to the conclusion and discussion. In the Appendix [A] we give 
the mass-squared matrices of the squarks to make our notation clear. 

II. GAUGE-HIGGS SYSTEM 

To discuss the statics and dynamics of the EWPT, we consider the gauge-Higgs system, 
which is governed by the Lagrangian 

/: g au g c-Hi gg s = -\FZ v F a » v - -B^BV + {D^ifD^t - V eS ($ d , $ u ; T). (2) 

i=d,u 

The covariant derivatives are 

D^ d = (dp + itfr^Al - *| V) *d, D^u = (d, + ig^Al + z| V) ^ (3) 

where r a is the Pauli matrices, #2 and g\ are the gauge couplings for SU{2)l and U(1)y, 
respectively. V e s is the effective potential, which is composed of the tree and one-loop parts 

V eS ($ d , $ u ; T) = V ($ d , $ u ) + AV($ d , $ u ; T). (4) 
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The tree-level part is given by the F-, D-, and the soft supersymmetry (SUSY) breaking 
terms 



U ($ d , $J = m\& d <S> d + m^t$ u - (mleij^i + h.c.) 



+ g -^(^ d ~ ^u) 2 + f (<fe)(<^), (5) 



where 



2 ~ 2 , | |2 2 ~ 2 , I |2 2 d /c\ 

m 1 = m 1 + \n\ , m 2 = m 2 + , m 3 = ^±5. (b) 
The Higgs doublets are parameterized as 

where we assume that U(l) em is not broken at the vacuum. Unlike the general two-Higgs- 
doublet model [la ]. CP violation cannot be accommodated in the Higgs potential at the 
tree level, leading to $ = 0. 

The one-loop part comprises the effective potentials at zero and nonzero temperatures 



T 4 . (m 2 



Fo(m\) + — 2 I B , F iji 



(8) 



AV($ dl $ u ;T)=^c A 

A 

where 



where Fq is the zero temperature effective potential, which is regularized in the DR scheme, 
ca denotes the degrees of freedom of particle species A, rfiA is the background- field-dependent 
mass, and M is a renormalization scale, which is determined in such a way that the zero 
temperature one-loop correction to the effective potential vanishes in the vacuum. In what 
follows, the contributions of the weak gauge bosons (Z,W), the third-generation quarks 
(t,b), and squarks (£1,2,^1,2) ar e taken into account. The statistical factors of each particle 
is, respectively, given by 

c z = 3, c w = 6, c t = c b = -AN C , Cfl2 = c- bi2 = 2iV c , (10) 

where N c is the number of color. In this article, we use the one-loop mass formulae of the 



Higgs bosons presented in Ref. 



a. 



III. LIGHT HIGGS BOSON SCENARIO 



To make a discussion on the LHS clear, we rotate the two Higgs doublets (ir 2 § d , $ u ) into 
(zr 2 $^*,$' u ) by j3 = tan~ 1 (v u /v d ). For simplicity, we assume that CP is conserved for the 
moment. The rotated Higgs doublets are cast into the form 




where t>o = yv \ + f 2 ~ 246 GeV, A is the CP-odd Higgs boson, H + is the charged Higgs 
boson and (G°, G ) are the Nambu-Goldstone bosons. h' d and h' u are related to the physical 
CP-even Higgs bosons (h, H) via 

/ 




(12) 



h' d \ I cos((3 — a) sin(/? — a) 
\h' u J \ — sin(/3 — a) cos(/3 — a) 

where a is a mixing angle between h d and h u . It should be noticed that only h' d possesses 
the Higgs VEV in the rotated basis, and we can identify it the SM-like Higgs boson, i.e., 
h' d = h SM . 

At the tree level, the Higgs boson couplings to Z boson normalized to the corresponding 
SM values are, respectively, given by 

9hZZ .fr, \ 9HZZ (a \ /io\ 

sm(p — a), = cos(p — a), (16) 



dh^zz 9hs^zz 

9hZA /a x 9HZA ■ , \ /-. A \ 

cos(p — a), = sm(p — a). (14) 



9h SM ZG° 9hSM Z G° 

From Eqs. f|TT]) . fTT2l . ffT3|) and ffl4l . one can see that h — > h SM , ghzz 9h SM zz an d 9hza 
for sin(/3 — a) — > 1 and _P — * /i SM , — > 9h sm zz an d fi'/fZA — * for cos(/3 — a) — > 1. 
The former case ("decoupling limit") is realized by taking rriz vtia and the latter one 
( "antidecoupling limit") is possible for rriA ~ which corresponds to the LHS. Depending 
on the theory parameters, some of the above four Higgs boson couplings can be small enough 
to avoid the LEP exclusion limits. In the LHS, for example, h can be as light as about 100 
GeV for ghzz/ 9h sm zz — 0.5. However, from the sum rules among the Higgs boson couplings, 
i- e -, (9hzz + 9hzz)/91sm zz = (glzA + 9hza)/91sm zg o = 1, not all Higgs bosons are allowed 
to be light. If the Higgs boson is the SM-like, the mass of it should be larger than 114.4 GeV 
since its coupling to the Z boson becomes the SM value. We will discuss the experimental 
constraints in more detail in Sec. IVII1 
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Finally, we comment on the CP- violating case. Although CP is conserved in the tree- 
level Higgs sector, it can be broken by the radiative corrections such as the top /stop loops. 
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20j. In this 



The realization of the maximal CP violation is called the CPX scenario 
scenario, it is possible to accommodate an even lighter Higgs boson. For instance, the mass 
of the lightest Higgs boson can be as small as 40 GeV, which is experimentally viable. It 
is interesting to discuss the phenomenological consequences of this scenario. However, the 
CPX scenario requires the relatively large ratios of (fi/rriq, A t /rriq) assuming rriq = mf , 
which is not favored in the context of the strong first order EWPT that we will discuss in 
the next section. 

IV. ELECTROWEAK PHASE TRANSITION 

To explain how the first order EWPT is strengthened in the MSSM, we here give a brief 
review on the light stop scenario. To simplify the discussion, we use the high temperature 



expansion in the effective potential 21j. In the decoupling limit, h is responsible for the 
electroweak symmetry breaking since h = h SM . On the other hand, v is mostly shared by H 
in the LHS. 

If the EWPT is first order, the critical temperature (Tc) is given by the temperature 
at which the effective potential has two degenerate minima. Let Vc be the nonzero VEV 
at Tc- From the argument using the high temperature expansion: a < 1 in Ib,f(^ 2 )-, one 
can see that vq/Tq can be enhanced when the quartic term in V e s gets smaller and/or the 
coefficient of the cubic term with a negative coefficient (denoted E) becomes larger. The 
former implies that the light h/H is favored in the decoupling/LHS limit. The latter can 
be realized by the additional contributions from the bosons. In fact, the effect of the lighter 
stop t\ on E is sizable due to the large top Yukawa coupling constant (y t ) and the degrees 
of freedom, i.e., 2N C . In the high temperature expansion, one can obtain 

V eS (v; T) 3 —ETv 3 ~ -(E SM + E h )Tv\ (15) 

where 

l/i ^ N c \y t \ 3 sin 3 (3 ( |X t | 2 \ 3/2 , x 

E SM ~—^ (2m 3 w + m 3 z -0.01, E h ~ cl ™ ~ 1 ~ Mr > 16 
47rfg 12V27T V mj J 

where X t = A* t —ji/ tan (3 and the 0(g 2 ) contributions are neglected in the stop contribution, 

also the hierarchies between the soft SUSY breaking masses are implicitly assumed, namely, 



trig 3> Xt, m i R = are taken to be consistent with the LEP bounds on the Higgs boson mass 
and the p parameter. As a result of mi = 0, which is necessary for the enhancement of the 
loop effect, the lighter stop mass becomes less than the top mass m^ < m t . From Eq. 0161) . 
one can easily see that the no- mixing case X t = maximizes the loop effect. If we assume 
tan/3 = tan/3(T = 0), we obtain Ef ~ 0.056, which is about 6 times as large as the SM 
contributions. 

To make the analysis on the EWPT precisely, the corrections from the daisy diagrams 
must be taken into account, which yield the temperature dependent terms in the field- 
dependent masses of particles. For the light stop, we effectively have the following mass 
shift to leading order: 

m\ R (T)=m} R + cT\ (17) 

where c is positive and composed of the relevant couplings in the theory. Therefore, the 
EWPT would get weaker even if m~ R = 0. Nonetheless, as advocated in Ref. 6j the choice 
of m-^ = — cT 2 , which can induce the charge-color-breaking (CCB) vacuum, may still realize 
the same enhancement of the stop loop effect as outlined above. For simplicity, we will not 
pursue this possibility and take m~ R ~ putting the daisy diagrams aside in this article. 
Our analysis on the strength of the EWPT is expected to be the same, or might be more 
optimistic than the more realistic case depending on the model parameters. 

The high temperature expansion makes it easy to see how the first order phase transition 
is strengthened analytically. It is, however, untrustworthy to use the approximation when 
the masses of the particles in loops are somewhat larger than Tq. We thus adopt a different 
method. Although the numerical integrations contained in the definitions of Ib,f{& 2 ) are the 
standard way, it is an extremely time-consuming task. Therefore, we will use the following 
fitting functions in our numerical analysis: 

N 

W(a 2 )=e- a E^ (18) 

n=0 

where the coefficients c b ^ are determined by the least squared method. For iV = 40, 
the errors of /^(a 2 ) do not exceed 10~ 6 for any a, which is sufficient for our purpose. 
To determine Tq and vo, we minimize V c s with Eq. (fl8l) in the three dimensional space 
(vd, v u cos v u sin •&) numerically. 
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V. CRITICAL BUBBLES 



The first order EWPT begins at somewhat below Tc- For the EWPT to proceed, the 
radius of the bubble should be larger than some critical size, otherwise it would shrink by 
the surface tension of the bubble wall. The bubble of this critical size is called the critical 
bubble. After nucleation of the critical bubbles, they start to percolate and eventually 
convert the symmetric phase into the broken phase if the supercooling is not too large. 

From Eq. (j2J), the energy functional in the temporal gauge takes the form 



E 



d 3 x 



1 /•""/••" + _R..R.. 



(19) 



(A$d) f A$d + (A$«) f A$« + Kff($ d , T) 

Here, we assume that the least energy has the pure-gauge configurations for A a (x) and B(x), 
hence = Bij = 0. Since a spherically symmetric configuration can give the least energy, 
the Higgs fields depend only on radial coordinate r = \fx?. The classical Higgs fields are 
parameterized as 

,ie d (r) ( j( r ) \ P i(8u(r)+p) / Q 

Pu{r) 



Pd[r) 




$«(r) 



(20) 



y/2 \ T ' v2 

where the gauge-invariant constant phase p will be taken to make the boundary value 9 a + 9 U 

at r = oo to vanish. To remove the energy source for the Z boson from the Higgs current, 

the so-called "sourcelessness condition" [22] 

2 d6d 2 d9u n 
P '-fr- p «-fr =0 (21) 

must be satisfied. Now let us introduce 9 = 9d + 9 U and 9 — 9d — 9 U . One can easily see that 
only 9 is a gauge-independent phase. Equation (12T1) can be rewritten as 



iP 2 d + Pl) Tr +(p 2 d -pl) Tr 



0. 



which is used to eliminate 9. The energy functional is then reduced to 



E 



Arc / dr r 2 
Jo 



dpd 
dr 



+ 



dpu" 2 
dr 



pIpI 



d9 



Y +V eS ( Pd ,p u ,9;T) 



2 P 2 d + Pl\drJ 

From this, the equations of motion (EOM) for p,i, p u and 9 are, respectively, given by 



, r 2^Pd 

r 2 dr\ dr 1 " 



LA 

r 2 dr 
1 d 



r 2 dr 



>dpu 
dr 



+ Pu 

1 d 

r dr 



Pi d&\ 


2 

+ 


dV cS 


Pd + Pl dr J 


dpd 


Pd d9\ 


2 

+ 


dV eS 


Pd + P 2 u dr) 


dp u 


2 PdPl dd 


) + 


dV eS 


Pd + Pi dr 


09 



0. 



(22) 

(23) 

(24) 
(25) 
(26) 



The solutions can exist only when the temperature lies in the range To < T < Tc, where 
To is the temperature at which the Higgs potential at the ori gin is destabilized in some 
direction. Since we focus on the CP-conserving bubble walls 23], 8(r) is set to be zero 
in the following discussion. The studies of the CP-violating bubble walls can be found in 
Refs. 
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241 ] . The boundary conditions for EOM are imposed in the symmetric (r = oo) 
and broken (r = 0) phases as 



lim p d (r) = 0, lim p u (r) = 0, 



dp d (r) 



dr 



0. 



dp u {r) 



r=0 



dr 



0. 



r=0 



(27) 



It is convenient to parameterize the Higgs profiles (p^, p u ) in terms of the dimensionless 
quantities 



Pd(r) 
v cos (3 ' 



Then E takes the form 



E = 4ttv 



dt e 



1 r /dht\ 2 
2\VdJ 



cos^ + (^ 



MO 



dh-: 



Pu(r) 
v sin p 



sm 2 p\ + V eH (h u h 2 ;T) 



(28) 



(29) 



where V e g = V e s/v 4 . Correspondingly, the EOM are rewritten as 



1 d f^dhx\ 1 



dV e 



cff 



£ 2 d£ V d£ J cos 2 (3 dhi 



1 ± ( 



,dho 



+ 



1 dV c 



cff 



£ 2 d£ V d£ J sin 2 (3 dh 



0. 
0. 



with the boundary conditions 



lim ^(0 = 0, Um/i 2 (0=0, 



dh^) 



d£ 



0. 



dh 2 {i) 



di 



o. 



5=0 



(30) 
(31) 

(32) 



We will comment on the numerical method to solve the EOM for /ii(£) and ^(0 i n Sec. 

EE 



A. Bubble nucleation temperature 



The bubble nucleation rate per unit time per unit volume is given by 25] 

3/2 

1 



T N (T) ~ T l (^f 1 ) e- E -< T » T . (33) 



where E ch (T) is the energy of the critical bubble at temperature T. Note that this is a rate 
per unit volume. We define the nucleation temperature T N as the temperature at which the 
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rate of the nucleation of a critical bubble within a horizon volume is equal to the Hubble 
parameter at that temperature. Since the horizon scale is roughly given by if _1 (T), the 
nucleation temperature is defined by 

T N (T N )/H 3 (T N ) = H(T N ) ~ 1.66Jg*(T N )T*/m P , (34) 



where g*(T N ) is the massless degrees of freedom at T N and m-p is the Planck mass (~ 
1.22 x 10 19 GeV). Since a single bubble nucleated within the horizon volume cannot convert 
the whole region to the broken phase, the nucleation temperature defined by Eq. (134"]) simply 
gives an upper bound of the temperature at which the EWPT starts. To determine T/v more 
precisely, simulations or the methods employed in Ref. are needed. 
From Eq. (|34|) . it follows that 

E ch (T N ) 3 ln /^v)\ =15259 _ 21n ^ (7V) _ 41n (^A_\ (35) 



T N 2 \ T N J " v ' V100 GeV 

Roughly, E c \y/T < 150 is necessary for the development of the EWPT. 

B. Numerical algorithm for critical bubbles 

To solve the EOM (1301) and (!3T!) under the boundary conditions (132]) . we use "relaxation 
methods." To implement the relaxation methods successfully, the initial configurations of 
the bubbles are of great importance. Here, we take the kink- type ansatz for them, i.e., 

kl ® = h2 ^ = l-tanh(-i^) ' (36) 
where R is the radius of the bubbles and L w is the bubble wall width. Since the best values of 

R and L w can vary as the temperature goes down, it is impossible to know them in advance. 

Nonetheless, it turns out that in most cases 10 < L w < R < 100 can give the successful 

initial configurations as long as T lies somewhere between T/v and Tc- After obtaining the 

convergent solutions at some temperature, we use them as the initial configurations to search 

for the critical bubbles at the temperature which is near the previous one. By doing this 

process iteratively, we find T N . 

VI. SPHALERON DECOUPLING CONDITION 



As shown by Manton in Ref. 27|, the static and unstable classical solution with a finite 
energy can exist in the SU (2) gauge theory. Such a configuration is called sphaleron, which 
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has a maximal energy along the least energy path. Similar to the instanton process at zero 
temperature, the baryon number can be violated by the sphaleron process at finite temper- 
ature. Although the instanton process, which is quantum tunneling, is highly suppressed 
and unobservable, the sphaleron process, which is thermal fluctuation, can be active at high 
temperature, especially before the electroweak symmetry is broken. 

In order to avoid the washout of the generated BAU after the phase transition, the 
sphaleron process must be decoupled. This condition can be obtained by demanding that 
the baryon number changing rate be smaller than the Hubble parameter 28 j 

1 dB 13N f uj- 



B dt 



4 • 32vr 2 a w 



KjV tI Kote- E ^ /T < H(T), 



(37) 



where Nf is number of generation, aw = ot em / sin 2 6{mzL oj- is the negative mode of the 
fluctuations around sphaleron, k is the 0(1) coefficient [3, 0]> -Mr an d Afct are contribu- 
tions from the translational and rotational zero modes, respectively, _E sp h is the sphaleron 
energy. If we denote the sphaleron energy as E sph = 4tcv£ / g 2 , Eq. (|3"7|) can be translated 
into 



T 4?r£ 



42.97 + lnfftMrJVrot) + In 



(-) 

\m w / 



( 9* 



-In . 

2 V 106.75 



-2 In 



T 



100 GeV 



.(38) 



It should be noted that model dependent parameters other than the sphaleron energy in the 
right-hand side of Eq. (|38|) contribute only logarithmically. 



£ 



In the SM 
= 2.00 



311 ]. the sphaleron so 



33j, MXot = 80.13 



29 



ution obtained with the tree-level Higgs potential yields 
24|, uj 2 _ = 2.3m w [J, |3J, ba] for \jg\ = 1. Inserting 



these values into Eq. ( 1381) together with k = 1, g* — 106.75, and T = 100 GeV, we find 



T 



> 0.026 x (42.97 + 4.38 + 0.416) = 1.24, 



(39) 



where g 2 = 2^Jira cm (mz)/ smOwijnz) = 0.652 ]lj] is used. It is found that the contributions 
of the zero mode factors to v/T can reach around 10%, while the contributions from the last 
three terms in the right-hand side of Eq. (1381) are only about 1%. 

Now we discuss the sphaleron decoupling condition in the MSSM. As an example, we 



take tan/? = 10.11, m H ± = 127 A GeV, A t = A b 



-300 GeV, n = 100 GeV, m. 



1200 



GeV, m t " = 10~ 4 GeV. This parameter set is experimentally allowed as we will discuss in 
Sec. IVIIII To see the effects of the temperature dependence of the sphaleron solution and 
the zero mode factors on the sphaleron decoupling condition, we consider the following three 
cases: 
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I: The sphaleron decoupling condition without the zero mode factors based on the zero 
temperature potential V e s{T = 0). 

II: The sphaleron decoupling condition with the zero mode factors based on the zero 
temperature potential V e s{T = 0). 

Ill: The sphaleron decoupling condition with the zero mode factors based on the full 
effective potential V e g(T ^ 0). 

The result of each case is summarized in Table I T = T/v is used in the right-hand side of 
Eq. (|38|) . As we see from I and II, the contribution of the zero mode factors to the sphaleron 
decoupling condition can be as large as 10% level, which is the same as in the SM case. 
In case III, vn/T n > 1.38 is obtained, which is about 40% stronger than the usual rough 
estimate used in the literature. We look into this case in more detail in the following. The 
temperature dependences of £, Af tT and Af rot are plotted in Fig. [TJ The endpoints of those 
curves correspond to those at Tq- The left panel shows that the sphaleron energy decreases 
as T increases. The change of the sphaleron energy is sizable around Tq, hence which can 
be the dominant error of the sphaleron decoupling condition if we neglect this effect. On the 
other hand, as shown in the right panel the temperature dependences of Af tT and J\f TOt are 
mild and thus have little effects on v/T. From this analysis, we conclude that in order to 
evaluate the sphaleron decoupling condition within a 10% accuracy, the sphaleron solution 
and the zero mode factors at finite temperature must be taken into account. 



Here, we comment on the earlier work 



361 ] . The authors have evaluated sphaleron at 



finite temperature using the one-loop finite-temperature effective potential with the daisy 
resummations. In their analysis, only the energy was calculated and the zero mode factors 
were missing. To search for the possible region where the EWPT is strong first order, they 



used the following sphaleron decoupling condition 371 ]: 

> 45. (40) 
In our calculation including the zero mode factors, this condition is translated into 

^ > 1-32 (41) 

in case III. 
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I 


II 


III 


£ 


1.89 


1.89 


1.77 


Mr 




7.36 


6.65 


Kot 




10.84 


12.27 


v N /T N > 


1.17 


1.29 


1.38 



TABLE I: The sphaleron decoupling conditions in the three cases I, II, and III. 
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FIG. 1: The temperature dependences of the sphaleron energy £ = E sp ^/ (At^v / g 2 ) and the zero 
mode factors Mr,rot- 

VII. EXPERIMENTAL CONSTRAINTS 

In this section, we present the experimental constraints that we will impose in the numer- 
ical calculation. When the masses of the three neutral Higgs boson (Hi) are smaller than 
114.4 GeV and/or the sum of two of them are smaller than about 195 GeV, we require 



9 2 Hi zz x Br(iZi -> //) < r Hi z( m H t ), 

9%h,z x Br(iZ< -> //) x Bx(Hj -> //) < TH iHj (m Hi + m H , 



(42) 



where f = b, r. ^HiZ, and J r H i H i are the 95% C.L. upper limits from the LEP experi- 



ments 
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14J . The Higgs boson couplings to the Z boson normalized to their SM values 



are given by, respectively 

9HiZZ = On cos (3 + 2i sin ft, 



(43) 
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QBiHiZ = (Ou0 3j - OtjOsi) sin/? - (0 2i Zj - 2j 3i ) cos/3, (44) 

where O is a 3-by-3 orthogonal matrix, which diagonalizes the mass-squared matrix of the 
neutral Higgs bosons. 

According to the recent results of the direct search for the SM Higgs boson at the Fermlab 
Tevatron, the mass range of the SM Higgs boson, 160 GeV < m h sm < 170 GeV has been 



. Since the upper bound of the lightest Higgs boson mass in the 



39] , the above excluded mass range is not relevant for the current 



excluded at 95% C.L. 
MSSM is around 135 GeV 
investigation. 

We also consider the constraint on the p parameter, which is a measure of the custo- 
dial SU{2) symmetry breaking. Ap = A^p + A H p < 0.002 is imposed in our calculation, 
where Agp and AhP are the contributions of stop/sbottom and the physical Higgs bosons, 
respectively. It turns out that as long as we take mi R <C m^, Agp would not exceed the 
current upper bound even for -C mr due to the suppressed couplings in the loop. In 
most of the region under investigation, the mass difference between one of the neutral Higgs 
bosons and the charged Higgs bosons is small, which implies the custodial SU(2) symmetry 
approximately exists, we then have AhP — 0. 

The experimental lower bounds on the masses of the SUSY particles are also taken into 
account, especially for the lighter stop, chargino (xf), an d neutralino (x?)> which is the 
lightest supersymmetric particle in our case, we impose that > 95.7 GeV, m x ± > 94 
GeV and m x o > 46 GeV |l|. 

Currently, the enormous data of the B physics observables are available, and especially 
B u — > ru T , B — > Xg'-f and B s — > p + p~ are relevant for the LHS. The averaged experimental 



40|: 



values of those processes are reported by HFAG 

Br(B u - rz/ T ) exp = 1.41+°^ x 10" 4 , (45) 
Br(B -> V s7 ) exp = (3.52 ± 0.23 ± 0.09) x 10" 4 , (46) 
Bi(B s -> p + p~) cxp < 0.23 x 10- 7 . (47) 

As mentioned in the Introduction, the first two decay modes can give severe constraints for 
the light charged Higgs bosons and the last one could be important for the light neutral 
Higgs bosons. However, it is worth to noting that to suppress the p-parameter corrections, 
the mass difference of the charged Higgs boson and one of the neutral Higgs boson must 
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be small. Therefore, the above three B decay modes should be simultaneously taken into 
account. 

At the tree level, the ratio of the two branching ratios of Bi(B u — > rz/ T ) M gg M /s M is given 



by 



4l| 



Br(B u -> tv t )sm V m H± I 

From the latest data 1] \V uh \ = (3.95±0.35) x 1(T 3 , f B = 200±20 MeV, which are involved in 
Br(B u — > rz/ r )sM, it follows that r# = 1.28 ±0.52. We will take the 95% C.L. interval range 
of r H . As for Br(£> — > X s ^) and Br(£> s — > we use the public code CPsuperH2.0 [20] 

to find the allowed region. 



VIII. NUMERICAL ANALYSIS 
A. Light Higgs boson scenario 



Before going to show the results of the EWPT, the phenomenologically allowed regions 
are considered. The strongest constraint in the LHS mostly comes from B — > X s j. It is 
well known that the contribution of the charged Higgs boson loop can be cancelled by that 
of the chargino loop. To make this mechanism work, the signs of /i and A are important. 
In Fig. [21 the 2a allowed region for B — > X s ^ is shown in the plane. We take 

\A\ = \A t \ = \A b \, Arg(A t ) = Arg(A b ) = vr, Arg(^) = 0, tan/3 = 12, m H ± = 130 GeV, 
rriq = 1200 GeV, m iR = 10~ 4 GeV, m~ bR = 1000 GeV, M l = 100 GeV, M 2 = 500 GeV. The 
region between the two red lines is allowed. Since the signs of A and \x are different from each 
other, X t cannot be vanishing. It means that the maximal effect of the stop thermal loop 
on the strength of the EWPT cannot be realized in this case. Moreover, the requirement of 
100 GeV < < \A\ makes the situation worse. Here, the lower bound of comes from 
experimental bound of the lighter chargino. Since the small X t is favored from the strong 
first order EWPT point of view, we take /i = 100 GeV, A = —300 GeV in the following 
discussion. 

Figure [3] shows the combined results of the LEP and B physics experimental constraints, 
and we overlay the rrih = 90, 100 GeV contours. The region above the blue dotted line is 
excluded by the B u — > tu t mode, and below the red dashed-dotted line, B — > X s ^ exceeds 
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100 200 300 400 500 600 700 800 9001000 

\A\ [GeV] 

FIG. 2: The allowed region of B — > X s j in the |/^|-|^4| plane, tan/3 = 12 and m#± = 130 GeV are 
taken. 




m H ± [GeV] 



FIG. 3: The allowed region in the LHS. We take A t = A b = -300 GeV, /i = 100 GeV, m ? - = 1200 
GeV, mi R = 10~ 4 GeV, and other input parameters are presented in the text. 

the experimental la limit. After taking the LEP constraints into account, the region, which 
is surrounded by the green curves, is phenomenologically allowed. In this region, the mass 
of the lightest Higgs boson is mostly between 90 GeV and 100 GeV, which are shown as the 
two magenta dashed curves, where the left one stands for = 90 GeV. The mass of the 
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CP-odd Higgs boson is approximately given by tua — rrih + 3 GeV and that of the heavy 
CP-even Higgs boson is slightly above the LEP exclusion limit, i.e., 114.4 GeV. For this 
parameter set, the constraint of B s — > p, + pT does not play a role unless tan/3 > 20. It is 
found that vc/Tq in the allowed region is larger than 0.9, and the maximal value of it is 
found to be 

vc = 107-10 GeV 

T c 116.27 GeV V ; 

for tan/5 = 10.1 and = 127.4 GeV. This value is far below the required sphaleron 

decoupling condition as shown in Table [H 

As we discussed in Sec. |V] the temperature at which the EWPT occurs is not Tc but 
somewhat below it. With the same input parameter set as is used in Eq. (1491) . we find the 
critical bubble solutions and evaluate their nucleation temperature. In Fig. HJ P c b(T)/P 
is shown as a function of T. The blue dotted line corresponds to E ch (T N ) /T N = 150.39, 
where Tn = 115.59 GeV. Since (Tc — T^)/Tc = 5.8 x 10 -3 , the supercooling is small. The 
left plot in Fig. [5] shows that the bubble walls at T = 116.00, 115.59 (= T N ), 114.00 GeV. 
The red solid curves represent foi(£), and the blue dotted curves correspond to h 2 (^). As T 
decreases, the energy of the broken vacuum becomes lower, correspondingly, and the radius 
of the bubble becomes smaller as it should be. It is found that the wall width is rather 
thick, which is consistent with the previous studies {23! . 



24|. In the standard mechanism of 



electroweak baryogenesis in the thick wall regime, the variation of tan f3(r) along the line 
connecting the broken phase and symmetric phase is crucial to the amount of the net BAU. 
Here, we define 

A(3(r)=(3(r)-(3(T), (50) 

where tan/5(T) = v u (T)/v d (T) and tan/3(r) = p u (r)/p d (r) = tan/3(r)/i 2 (0//ii(0- A P( r ) 
at Tat is shown in the right plot in Fig. [51 This shows that the behaviors of pd{r) and p u (r) 
in the intermediate of r are slightly different. We observe that pd(r) approaches to zero 
somewhat faster than p u (r) does as r increases. Such an enhancement might play a role in 
generating the BAU. 

At T N , v N /T N is given by 

v N 116.73 GeV 



T N 115.59 GeV 



1.01. (51) 
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FIG. 4: E ch (T)/T as a function of Tin the LHS. The blue dotted line shows E cb (T N ) /T N = 150.386, 
where T/v = 115.59 GeV. We take tan/3 = 10.1 and m#± = 127.4 GeV and other parameters are 
the same as in the Fig. [3] 

One can see that there is the (9(10)% enhancement in v/T compared to that at Tq- However, 
as we discussed in Sec. EH the sphaleron decoupling condition at T N is given by 



Therefore, the sphaleron process is not decoupled at T/v either in this parameter region. We 
search for the maximal value of vc/Tq scanning tan/3 and for rriq = 1300, 1400 and 
1500 GeV. Our findings in the LHS are summarized in Table [III Every case is more or less 
the same as the first one discussed here. We thus conclude that there is no region where the 
sphaleron decoupling condition is satisfied in the LHS. 

B. The decoupling limit 

In the decoupling limit, h is the SM-like Higgs boson and must be heavier than 114.4 
GeV. The constraints from B physics are not so strong in this region. We then consider the 
no-mixing scenario, X t = and Af, = A t are taken. Other parameters are the same as those 
in the LHS discussed in the previous subsection. For rriq < 1600 GeV, there is no allowed 
region such that rrih > 114.4 GeV. As an example, we take irtq = 1700 GeV and scan tan/3 
and 77iff± in the ranges: 1 < tan/3 < 50, 100 GeV < < 1000 GeV. The maximal value 



^ > 1.38. 



(52) 
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£ r = [GeV" 1 ] 



FIG. 5: Left: hi(£) (red solid curve) and /12(C) (blue dotted curve) are plotted for T = 
116.00,115.59 (= Tat), 114.00 GeV. Right: Af3(r) as function of r = £/v N . The input parame- 
ters are the same as in Fig. [J] 



of Vc/T c is realized in the case of tan/3 = 42.6 and m H ± = 1000 GeV. v/Tat T c and T N 
are found to be 

v c _ 114.46 GeV v N _ 121.45 GeV _ 

7b " 116.99 GeV ~ °' 95 ' ~ 116.22 GeV ~ L ° 5 ' (53) 

Similar to the LHS, v/T is ameliorated at Tjv by about 10%. However, it is still not strong 
enough to satisfy the sphaleron decoupling condition which is given by /T N > 1.38. 

We search for the allowed region varying nig = 1800, 1900, 2000 GeV and investigate 
vc/Tc and vn/Tn- The results are summarized in Table HTT1 There is no significant change 
between the first case discussed above and the others. We also scan in the region, trig > 2000 
GeV. The larger rrig allows the smaller tan/3 region, which implies that the top Yukawa 
coupling becomes stronger. On the other hand, rrth also becomes larger through the one- 
loop top/stop corrections. The numerical results show that the first order EWPT is getting 
weaker and weaker as rriq increases. 
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IX. CONCLUSION AND DISCUSSIONS 



We have reanalyzed the strong first order EWPT in the MSSM considering the experi- 
mental results such as the LEP and the latest B physics data. v/T was evaluated at both 
T c and T N in the LHS and ordinary decoupling limit. In the LHS with a light stop, the 
no- mixing scenario, which can maximize the stop thermal loop effect, cannot be realized 
due to the severe constraint from B — > X s ^. Combining with the other experimental data, 
especially the LEP and B u — > tv t mode, the allowed region becomes more constrained. 

The sphaleron decoupling condition was also calculated using the one-loop effective po- 
tential at finite temperature. We found that the contributions of the zero mode factors 
coming from the fluctuations around sphaleron can be as large as about 10%. In the phe- 
nomenologically allowed region, the typical values of the sphaleron decoupling conditions at 
T/v are found to be v^/Tfq > 1.38. 

It is observed that about 10% enhancement of v/T is possible at T/v in comparison with 
that at Tc- We find that vn/Tn — 1.01 in the LHS and vn/Tn — 1.05 in the decoupling limit. 
However, such a strength of the EWPT is not enough to satisfy the sphaleron decoupling 
condition. 

Although electroweak baryogenesis in the MSSM seems to be infeasible, it may be possible 
to circumvent the above negative conclusion. Some comments on the possible loopholes are 
in order. Firstly, the sphaleron decoupling condition is supposed to be imposed at the 
temperature at which the EWPT ends, which can relax the decoupling condition. Secondly, 
our calculations, which are based on the one-loop effective potential, can be modified by the 
higher-order loop contributions such as the two- loop effects |42j. Lastly, as we mentioned 
briefly in Sec. [IV]the scalar potential can be extended in such a way that the lighter stop also 
has the VEV, which can develop the CCB vacuum. In Ref. it is claimed that the EWPT 
is strong first order when the electroweak vacuum is metastable and the CCB vacuum is 
the global minimum. Since the time scale of the decay of the metastable vacuum is longer 
than the age of the Universe, this scenario is viable. The two-loop contributions are also 
included in their calculation. It is interesting to investigate whether or not the EWPT is 
still strong enough to satisfy the sphaleron decoupling conditions that we have considered 



in this article. A devoted study of the EWPT and the spha' 



on the two- loop effective potential will be given elsewhere 43]. 



eron decoupling condition based 
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APPENDIX A: THE MASS-SQUARED MATRICES OF STOP AND SBOTTOM 



Here, we give the expressions of the mass-squared matrices of the stop and sbottom 

' m\ + + | (gl - f ) {vl - vl) %(A* t e-^v u - nv d ) 

v %A*«v u -fv d ) ml + ^fvl+ 9 4(vj-vl) 



(Al) 



^ = I " z u °r a ^ " - ^ " ~ ' | , (A2) 

where y t and are the top and bottom Yukawa couplings, respectively. 
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TABLE II: Examples of the first order EWPT in the LHS. \A b \ = \A t \ = 300 GeV, m~ tR = 10~ 4 
GeV, m~ bR = 1000 GeV, = 100 GeV, M x = 100 GeV, M 2 = 500 GeV, Arg(A t ) = Arg(A 6 ) = 
7T, Arg(/i) = 0. 



niq (GeV) 


1200 


1300 


1400 


1500 


m h (GeV) 


164.76 


165.66 


166.37 


166.95 


m ± (GeV) 


94.86 


94.80 


94.77 


94.72 


m x o (GeV) 


65.14 


65.04 


64.99 


64.92 


m x o (GeV) 


108.31 


108.26 


108.24 


108.21 


m x o (GeV) 


129.46 


129.50 


129.52 


129.56 


tan (3 


10.11 


9.87 


9.75 


9.57 


m H ± (GeV) 


127.40 


127.40 


127.50 


127.50 


m Hl (GeV) 


94.04 


93.95 


94.03 


93.97 




Q7 89 
y 1 .oZ 


Q7 


Q7 QR 

y i .yo 


Q7 QQ 

y I .yy 


m# 3 (GeV) 


116.47 


117.13 


117.72 


118.31 


9h-\VV 


0.228 


0.228 


0.227 


0.227 


9h 2 vv 


0.000 


0.000 


0.000 


0.000 


9h 3 vv 


0.772 


0.772 


0.773 


0.773 


vc/Tc 


107 096 

= 0.921 

116.274 


107 512 

= 0.923 

116.496 


107 769 

: = 0.923 

116.770 


107 915 

: = 0.922 

117.045 


tan /3c* 


13.803 


13.640 


13.597 


13.455 




1 1 fi 797 

Lib - 727 = 1.010 
115.585 


1 1 7 1 ^ 

ii7 - i55 = 1.012 
115.798 


117 4D4 
117.404 = 

116.067 


1 1 7 1 

Ln - b6L = 1.010 
116.339 


tan /?tv 


13.676 


13.503 


13.453 


13.307 


E c b(T N )/T N 


150.386 


150.379 


150.370 


150.360 


£ 


1.769 


1.770 


1.770 


1.771 


AA tr 


6.652 


6.658 


6.662 


6.667 


A/" rot 


12.266 


12.253 


12.240 


12.229 


W^/v > 


1.383 


1.382 


1.382 


1.380 
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TABLE III: Examples of the first order EWPT in the no-mixing scenario. \At\ = \Ab\ = \fi\/ tan (3, 
m~ tR = 10" 4 GeV, m~ bR = 1000 GeV, = 100 GeV, Mi = 100 GeV, M 2 = 500 GeV, Arg(A t ) = 
Avg(A b ) = Arg(^) = 0. 



mq (GeV) 



m h (GeV) 
to + (GeV) 
to x o (GeV) 
m x o (GeV) 
to x o (GeV) 



tan/? 
m H ± (GeV) 
m Hl (GeV) 
TOff 2 (GeV) 
m# 3 (GeV) 

SHiVV 
9h 2 VV 
9h 3 VV 



vc/T c 
tan (3c 
vn/T n 
tan /3tv 
E ch (T N )/T N 



£ 
AA tr 

A/" rot 

t>Ar/T/v > 



1700 



170.73 
96.80 
68.25 
109.73 
128.15 



42.62 
1000.00 
114.40 
994.09 
994.11 
1.000 
0.000 
0.000 



111.461 



116.993 

42.966 
121.454 



0.953 



116.221 
42.955 

150.366 



1.045 



1.773 
6.677 
12.211 
1.379 



1800 



170.76 

95.69 

66.46 

108.92 

128.91 



15.10 
1000.00 
114.42 
996.45 
996.54 
1.000 
0.000 
0.000 



111.460 



117.007 

15.171 
121.452 



0.953 



116.236 
15.168 

150.370 



1.045 



1.773 
6.677 
12.210 
1.379 



1900 



170.78 

95.06 

65.45 

108.45 

129.34 



10.97 
1000.00 
114.40 
996.61 
996.77 
1.000 
0.000 
0.000 



111.483 



116.994 

11.022 
121.478 



0.953 



116.222 
11.019 

150.364 



1.045 



1.773 
6.678 
12.210 
1.379 



2000 



170.81 
94.66 
64.82 
108.16 
129.59 



9.35 
1000.00 
114.48 
996.66 
996.88 
1.000 
0.000 
0.000 



111.440 



117.060 

9.394 
121.424 



0.952 



= 1.044 



116.288 
9.392 

150.360 



1.773 
6.678 
12.209 
1.379 
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